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Abstract 

l/") . In this note two cosmological models representing the flat Friedmann Universe filled with a 

0^ ■ Chaplygin fluid, with or without dust, are analyzed in terms of the recently proposed "statefinder" 

parameters [ffl. Trajectories of both models in the parameter plane are shown to be significantly 
different w.r.t. "quiessence" and "tracker" models. The generalized Chaplygin gas model with an 
equation of state of the form p = ~A/p a is also analyzed in terms of the statefinder parameters. 
PACS numbers: 98.80.Es, 98.80.Cq, 98.80. Hw 
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In the search for cosmological models describing the observed cosmic acceleration 
inspiration coming from inflation has suggested mainly models making use of scalar fields 
There are of course alternatives; in particular, in JToj , [Tl| , p"2| an elementary model has been presented 
' describing a Friedmann universe filled with a perfect fluid obeying the Chaplygin equation of state 
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p-$ (1) 



where A is a positive constant (for a thorough review see Ref. (T^j). The interesting feature of this 
model is that it naturally provides a universe that undergoes a transition from a decelerating phase, 
driven by dust-like matter, to a cosmic acceleration at later stages of its evolution (see for details). 
An interesting attempt to justify this model H| makes use of an effective field theory for a three-brane 
universe [fT5| . 

In the flat case, the model can be equivalently described in terms of a homogeneous minimally 
coupled scalar field 0, with potential JTc| ] 

V{<j)) = \ VA (cosh 30+ — — I. (2) 



2 \ cosh 30 

However, since models trying to provide a description (if not an explanation) of the cosmic acceleration 
are proliferating, there exists the problem of discriminating between the various contenders. To this 
aim a new proposal introduced in Jlj makes use of a pair of parameters {r, s}, called "statefinder" . The 
relevant definition is as follows: 

r ^^' S = 3(<?-l/2)' (3) 

where H = - is the Hubble constant and q = —-^jp is the deceleration parameter. The new feature of 
the statefinder is that it involves the third derivative of the cosmological radius. 
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Trajectories in the {s, r}-plane corresponding to different cosmological models exhibit qualitatively 
different behaviors. ACDM model diagrams correspond to the fixed point s = 0, r = 1. The so-called 
"quiessence" models Q are described by vertical segments with r decreasing from r = 1 down to some 
definite value. Tracker models ]f8| have typical trajectories similar to arcs of parabola lying in the 
positive quadrant with positive second derivative. 

The current location of the parameters s and r in these diagrams can be calculated in models 
(given the deceleration parameter); it may also be extracted from data coming from SNAP (SuperNovae 
Acceleration Probe)-type experiments B. Therefore, the statefinder diagnostic combined with future 
SNAP observations may possibly be used to discriminate between different dark energy models. 

In this short paper we apply that diagnostic to Chaplygin cosmological models (the direct comparison 
of the already available supernova data with the Chaplygin gas model was also undertaken recently in 
[jig , [l7|]). We consider both the one- fluid pure Chaplygin gas model and a two- fluid model where dust 
is present as well. We show that these models are different from those considered in |[| and they are 
worth of further study. 

To begin with, let us rewrite the formulae for the statefinder parameters in a form convenient for 
our purposes. We shall need the Friedmann equation for the flat universe 



H 2 = ~ 2 =P (4) 



and the energy conservation equation 



p=-3H(p + p). (5) 
Using these two equations it is easy to find that 

2 2 p 

and then 

r = 1 - _?E_ s = t- ( 7 ) 

For a one-component fluidp] these formulae become especially simple. Since 

dp . dp 
- p = -3 V p {p+p) - 



P=^P=-3VP(P+P)^, (8) 



we easily get: 

'- 1+ §K)£- »=H)t- <9) 

For the Chaplygin gas one has simply that 

v>= d / = ^ = - P - = l + s (10) 
op p A p 

(v^ is the square of the velocity of sound) and therefore 

r = l-|*(l + *). (11) 

Thus, the curve of r(s) is an arc of parabola. To find the admissible values of s, we note that Eqs. ([!]), 
(pj) and (||) easily give the following dependence of the energy density on the cosmological scale factor 



(12) 



1 We confine ourselves to the case, which is fulfilled in our models, of a fluid for which the equation of state has the 
form p = p(p). 



3 



where B is an integration constant; therefore 



«; = i+* = ^-t-b- ( 13 ) 
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When the cosmological scale factor a changes from to oo the velocity of sound varies from to 1 and 
s varies from —1 to 0. Thus in our model the statefinder s takes negative values; this feature is not 
shared by quiessence and tracker models considered in |j| . 

As s varies in the interval [—1,0], r first increases from r = 1 to its maximum value and then 
decreases to the ACDM fixed point s = 0, r = 1 (see Fig. pi). 




Figure 1: s-r evolution diagram for the pure Chaplygin gas 



If q « —0.5 the current values of the statefinder (within our model) are s w -0.3, r w 1.9. In § it is 
reported an interesting numerical experiment based on 1000 realizations of a SNAP-type experiment, 
probing a fiducial ACDM model. Our values of the statefinder lie outside the three-sigma confidence 
region displayed in Based on this fact it can be expected that future SNAP experiments should be 
able to discriminate between the pure Chaplygin gas model and the standard ACDM model. 

Let us consider now a more "realistic" cosmological model which, besides a Chaplygin's component, 
contains also a dust component. For a two-component fluid Eqs. @) take the following form: 



2(p + pi) 



dp 
dp 



(P + P) 
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dpi 
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1 \dp, dpx 

s =— fl-(P+P) + ^— \Pi +Pi) 

p + pi [op dpi 

If one of the fluids is dust, i.e. pi = Pd = 0,the above formulae become 

_ p + p dp 
p dp 
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9(p + p) dp 
2(p + p d ) dp 



(14) 



(15) 



(16) 
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If the second fluid is the Chaplygin gas, proceeding exactly as before we obtain the following relation: 



r = l- 9 - S ±±2l (17) 



p 

To find the term pn/ p we write down the dependence of the dust density on the cosmological scale factor: 

Pd = S, (18) 
where C is a positive constant. Eq. ([l3| ) gives that Aa 6 + B = — — and therefore 

Pd C 



p y'Aa 6 + B 



kV-s, (19) 



where the constant k — C/y/B is the ratio between the energy densities of dust and of the Chaplygin 
gas at the beginning of the cosmological evolution (cf. Eq. (03)). Thus 



Graphs of the function (20) for different choices of K are plotted in Fig. |2j. 

In this case there are choices of the parameters so that the current values of the statefinder are close 
to the ACDM fixed point; indeed for n — 1 we have s = —0.09 and r = 1.2835. By increasing k we 
get closer and closer. Already for k = 2 we get s = 0.035, r = 1.11 while for n > 10 the statefinder 
essentially coincides with the ACDM fixed point (see Fig. ||). 

Thus our two-fluid cosmological models (with n say bigger than 2) cannot be discriminated from 
the ACDM model on the basis of the statefinder analysis. On the other hand this fact could also 
be interpreted as an argument in favor of the model because the ACDM model does not contradict 
observations up to now. 

The attractive point of our two-fluid cosmological model with a parameter k of order one is that 
it may suggest a solution to the cosmic coincidence conundrum: indeed here the initial values of the 
energies of dust and of the Chaplygin gas are of the same order of magnitude. 

The comparison of the two-fluid (i.e the Chaplygin gas plus dust) cosmological model with observa- 
tional data has also been studied in |l6|, [l7| . An analysis of the data obtained from 26 Supernovae leads 
to the conclusion JT(| that the best fitting model seems to be the pure Chaplygin gas without dust. 
On the other hand |l7fl by combining the results of the observations of 92 Supernovae with the matter 
power spectrum a wide range of values of k seems to be admissible. 

Thus, for getting more precise constraints on the parameters of the Chaplygin models one needs 
both new observations and additional diagnostic techniques. From this point of view an application of 
the statefinder analysis could be useful. 

Similar conclusions also hold for a one-fluid model of a generalized Chaplygin gas with a modified 
equation of state, as introduced in p0[ : 

A , \ 

P = - — > 21 
P 

with < a < 1. This gives a cosmological evolution from an initial dust-like behavior to an asymptotic 
cosmological constant, with an intermediate epoch that can be seen as a mixture of a cosmological 
constant with a fluid obeying the state equation p = ap (a — corresponds to the ACDM model). This 
generalized model was studied in some detail in JTgj] . Equation ([H]) gives the following dependence of 
the energy density p on the scale factor a: 

B \ 
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Figure 2: s-r evolution diagram for the Chaplygin gas mixed with dust. Dots locate the current value 
of the statefindcr 



In this case the squared velocity of sound 

= - — (2>!) 



2 dp ap Aa 
dp P At a3{a+ i) 



varies from to a. From Eq. ([)]) it follows that s = v% — a and therefore the admissible values of s are 
now in the interval between —a and 0. From Eq. (j^) we get 

r = 1 _9_s( S Ta) 
2 a 

For small values of a the generalized Chaplygin gas model becomes indistinguishable from the standard 
ACDM cosmological model. 

In conclusion. From the statefinder viewpoint both the pure Chaplygin model and the two-fluid 
mixture have different behaviors w.r.t. other commonly studied models. In particular, a future larger 
amount of data on high z type la supernovae may allow to distinguish between the pure Chaplygin 
gas and the ACDM models. Instead, the mixed two-fluid model becomes practically indistinguishable 
from ACDM for sufficiently large values of the parameter k (k > kq ps 5). On the other hand a mixed 
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two-fluid model with k of order one is attractive from the point of view of a possible solution of the 
cosmic coincidence problem. 
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